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MICHAEL ANSHELEVICH 

Abstract. The operator L M : / i— > J £M^£il) d(j,(y) is, for a compactly supported measure ji 
with an L 3 density, a closed, densely defined operator on L 2 (n). We show that the operator Q = 
pL' 1 —qLp has polynomial eigenfunctions if and only if /j, is a free Meixner distribution. The only time 
Q has orthogonal polynomial eigenfunctions is if fi is a semicircular distribution. More generally, the 
only time the operator p(L„L M ) — qL^ has orthogonal polynomial eigenfunctions is when /i and v 
are related by a Jacobi shift. 
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1. Introduction 

The Meixner class of orthogonal polynomials, and its corresponding family of distributions — 
Gaussian, Poisson, gamma, binomial, negative binomial, and hyperbolic secant — has a number 
of combinatorial and probabilistic characterizations. For example, in the original paper [Mei34], 
Meixner showed that this class consists exactly of all orthogonal polynomials whose exponential 
generating functions have the special form A(z)e xB<yZ \ For probabilistic characterizations, see for 
example [LL60, Wes93]. This class also plays an important role in statistics, under the name of 
quadratic exponential families, see [Mor82, Mor83, DKSC08]. 

On the other hand, historically, the most important class of orthogonal polynomials, often called the 
"classical orthogonal polynomials" [AS90], consists of the Hermite, Laguerre, and Jacobi families. 
This class also has numerous characterizations, notably the following one in terms of second order 
differential equations [Ism05, Chapter 20]. 

Theorem (Bochner 1929). The Sturm-Liouville operator 

p(x)y" + q(x)y' 

has a family of eigenfunctions consisting of orthogonal polynomials if and only if these eigenfunc- 
tions are Hermite, Laguerre, or Jacobi. In this case p is a polynomial of degree at most 2, and q a 
polynomial of degree at most 1. 

In all of these cases, the operator is symmetric with respect to the inner product induced by the 
orthogonality measure of the polynomials. This orthogonality measure has a density determined by 

w' q — p' 



w p 

The class of probability densities with — = ''"f^. is larger than the Bochner class, and is known in 

r j w quadratic ° ' 

statistics as the Pearson class [DZ91]; the corresponding operator has polynomial eigenfunctions, 
which however are not necessarily orthogonal on the real line. See Remark 1 for more details. Note 
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that all measures in this class have continuous densities, while the Meixner class includes purely 
atomic measures such as the Poisson distribution. 

With motivation from free probability [VDN92, NS06], in [Ans03] I introduced the free Meixner 
distributions as all the probability measures whose orthogonal polynomials have the ordinary gener- 
ating function of the special form A(z) 1 _ rj } B ^ z y The free Meixner class also has an explicit descrip- 
tion, and is in fact quite simple, simpler than the classical Meixner class [CT84]. Nevertheless, the 
majority of characterizations of the classical Meixner class have their counterparts — some straight- 
forward, some not — for the free Meixner class [BB06, Bry09, Ans07, BL09]. The goal of this 
paper is to investigate the free version of the Bochner / Pearson class. We will show that, in the 
appropriate framework, the Bochner-type characterization holds for the free Meixner class itself, 
quite unlike in the classical situation. 

What is the "free version" of the derivative operator? Unfortunately (or, rather, fortunately), there is 
no classical-free theory dictionary. However, Voiculescu (in [Voi98, VoiOO] and subsequent work) 
has shown that some roles played by the derivative operator are played in free probability by the 
difference quotient operator, 

d : C[l] — > C[R x R], (df)(x,y) = Z^Ll/M. 

x-y 

It is hard to talk about differential equations in the difference quotient operator, since its domain 
and range are not the same. One modification of it is the single variable operator 

w = ^^>. 

x 

which is the q = version of the g-derivative operator from special function theory. The answer 
for the Bochner-type question for the g-derivative D q is known, see [AS90] or [Ism05, Chapter 
20]. However, for q — 0, all these classes of polynomials reduce simply to monomials (with the 
exception of the degree one polynomial), and so are quite uninteresting. Other questions of this 
type involve the Askey- Wilson operator V q , and the difference operator. In the latter case, the 
corresponding orthogonal polynomials are Charlier, Meixner and Krawtchouk (that is, the discrete 
measures in the Meixner class), but also the Hahn polynomials. Note that the Meixner-Pollaczek 
polynomials belong to the Meixner class without satisfying either a differential or a (real) difference 
equation. 

Instead of the choices above, we will consider the following single variable operator: for a measure 

H, 



L»[f] = {I®n)[df]= [ f(x) m d»(y). 

Jr x-y 



So this is really an integral operator, but we will think of it as an analog of the derivative opera- 
tor. Of course, the big difference from the derivative or the D operator is that it depends on \i. 
(Note, however, that D = Ls .) Our Bochner-type question is: for which ji can a second order 
operator in L M with polynomial coefficients have polynomial eigenfunctions? And when are these 
eigenfunctions orthogonal? 

Here is the summary of the main results of this paper. Section 2 contains, in particular, the detailed 
description of the Bochner and Pearson classes on one hand, and of the free Meixner class on the 
other, as a convenient reference. In Section 3, we show that an operator 

pLl + qL^ 
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has polynomial eigenfunctions if and only if /j, has a free Meixner distribution (given fi, the choice 
of p, q is unique). On the other hand, if p — 1 and q = —H, the conjugate variable for fi, then the 
operator pL 2 ^ + qL^ is symmetric. H is a polynomial of degree one exactly when \i is a semicircular 
distribution, which are also the only cases when the polynomial eigenfunctions are orthogonal with 
respect to /i. Finally, in Section 4 we look at operators of the form + qL^ which depend 

on a pair of measures. We classify pairs (/i, v) which correspond to polynomial eigenfunctions. 
For any /x, if v is a Bernoulli distribution, one can choose the coefficients in the operator so that it 
will have polynomial eigenfunctions. On the other hand, for any v, if fi is a "Jacobi shift" of v, the 
operator has polynomial eigenfunctions. Finally, we show that the eigenfunctions are orthogonal 
with respect to /i if and only if the pair (/i, v) is of the latter type; in that case, the operator is also 
symmetric. 

Acknowledgements. Thanks to Alberto Griinbaum for questions that motivated the writing of 
this paper, and to Serban Belinschi and David Damanik for useful conversations. Thanks also to 
Dennis Stanton, the person who knows everything about orthogonal polynomials, new and old, and 
generously shares it. 

2. PRELIMINARIES 

Notation 1. A polynomial system is a family of polynomials {P n : n > 0} with deg P n = n. 

We start with a more detailed description of Bochner's theorem, to be compared with the results 
that follow it. 

Remark 1 (Bochner- Pearson classification). If the operator 

y i-> (p(x)y')' + q(x)y' 

has a polynomial system of eigenfunctions, then p is a polynomial of degree at most 2, and q a 
polynomial of degree at most 1, and the operator is symmetric with respect to the inner product 
induced by the measure dfx(x) = w(x) dx with the density determined by 

w' q 
w p 

By affine transformations, one can reduce the analysis to the following cases. 

(a) ^ = —ax. Then a > (and we may taken a = 1), and the Hermite polynomials are 
eigenfunctions, orthogonal with respect to the Gaussian distribution -^e~ x I 2 . 

(b) ^ = ^ — (3. Then (5 > (and can be taken (5 = 1), a > 0, and the Laguerre polyno- 
mials are eigenfunctions, orthogonal with respect to the gamma distribution fjaj xa ~ le ~ x on 
(0,oo)._ 

(c) — = ^— ^ — jz^. Then a, (3 > 0, and the Jacobi polynomials are eigenfunctions, orthogonal 
with respect to the beta distribution I r ^p^ x Q! ~ 1 (l — x)^^ 1 on (0, 1). 

(d) ^ = -i±SE + a > o, > so take f3 = 1, and the measure is ^x'^e' 1 ^ on (0, oo) (Z 
will always denote the appropriate normalization constant). Bessel polynomials, orthogonal 
in the complex plane, are eigenfunctions. 

(e) ^ = + J^, a > 0, and the measure is \[x 2 + i)-(i+«)/2 e /Jtan-i(:0 For ^ = , 
these are the t-distributions. 
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Remark 2. While there is, in general, no relation between the Pearson and Meixner distributions, 
some parallels between these classes can be explained by the following observation. As mentioned 
above, for the Pearson class dfi(x) = w(x) dx, 

w'(x) d + ex 

w(x) 1 + bx + ex 2 

On the other hand, let {P n (x, t)} be the monic orthogonal polynomials for the measure ii (t is the 
convolution parameter), and denote by F their exponential generating function: 

F(x,t,z) = y2-P n (x,t)z n . 

H is in the Meixner class if and only if 



nl 

n=0 



d z F x — tz 



F ' 1 + bz + cz 2 ' 
Notation 2. For a probability measure /i, its Cauchy transform is 

Gfi(z) = / — - — dti(x). 
Jr z x 

Recall that /i can be recovered from it as the weak limit 

da(x) = lim Im G(x + is) dx. 

7T e^0+ 

The Hilbert transform of /i, defined Lebesgue almost everywhere, is the function 

H\u](x) = - lim Re G(x + is). 

7T e^0+ 

If /j, has a density in L P (R, dx),p> 1, then H\fi] e L P (R, dx) [Tri85, Chapter 4]. 
The i?-transform of \i is determined by 

G (r(z) + ^=z; 

it is an analytic function on a domain, but in this paper we will identify it with its power series 
expansion. 

Throughout most of the paper, /i (and later, also v) will be a compactly supported probability 
measure. In particular, its moments 

m n = n[x n ] 

are finite, and we will also consider /i as a linear functional on polynomials. In this case, its moment 
generating function is 

OO 

M(z) = M fl (z) = J2 m nZ n - 

n=0 

Remark 3. For a compactly supported probability measure fx with an L 3 (IR, dx) density, its renor- 
malized Hilbert transform 

H = H^ = 27r#[/i] 
is well defined and is in I/ 2 (/x). Moreover, for / e L 2 (/j,), 

(jM®H)[df\=IM[Hf\. 
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In other words, H is the conjugate variable for \x; see [Voi98] for more details. As explained there, 
the conjugate variable is the free analog of the classical score function — ^ of a probability measure 

w(x) dx. 



Remark 4 (Free Meixner distributions). The normalized free Meixner distributions fj,b,c are proba- 
bility measures with Jacobi parameter sequences 

{(0,6,6,6,...),(l,l + c,l + c,l + c,...)}, 

for b E R, c > —1. The general free Meixner distributions are affine transformations of these. More 
explicitly, the distribution with parameters b, c is 



i ^( 4 (1 + c)- (* -&)*)_ 



dx + 0, 1, or 2 atoms. 



2tt 1 + bx + ex 2 

Unfortunately, none of the descriptions of these measures in [CT84, SY01, Bry09], including our 
own [Ans03], are complete, so we provide a detailed description here. By affine transformations, 
the situation can be reduced to the following six cases. For future reference, we also record their 
conjugate variables. 

i. b = c = 0, 

dfx(x) = — a/(4 — x 2 ) + dx. 
2n 

This is the semicircular distribution. Conjugate variable x. 

ii. c = 0, b j£ 0, a > 0, 



! J((l + V^) 2 -x) (x-(l-v^) 2 ) + 

dx + max(l — a, 0)5o. 

2n x 

These are the Marchenko-Pastur distributions. Conjugate variable — - + 1. 
iii. -1 <c<0, a,P>0, a + P> 1, 



i */(4a(l - x) - (a + 1 - (a + /3):r) 2 ) 



<ix + max(l — a, 0)5 + max(l — f3, 0)5i. 



2ir x(l — x) 

These are the free binomial distributions, including the Bernoulli distributions for c = — 1, 
a + (3 — 1. Conjugate variable — • 
iv. c ^ 0, 6 2 - 4c = 0, a > 0, 



1 ^((1 + V^TT) 2 - axj (ai - (1 - v/cTTT) 2 ) 

2^ x 2 
Conjugate variable 2j± ^ — 
v. c> 0, b 2 - 4c < o/a >0,(3eR, 



dx. 



J (a 2 x 2 - 2(3(2 + a)x + 4(/3 2 - 1 - a)) 

— dx. 



2tt 1 + x 2 

Conjugate variable - 
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vi. c> 0, b 2 - 4c> 0, a + (3 < 0, a/3 < 0, 



y (4a(l - x) - (a + 1 - {a + f3)x) 2 ^ 



dx 



2ir x(x — 1) 

+ ((1 - a)S for < a < l) + ((1 - for < f3 < l) ; 

note that since a + /3 < 0, at most one atom may appear. Conjugate variable ^f- — ^5§. 

Using the (6, c) parametrization, and either the orthogonal polynomials [Ans07] or the Laha-Lukacs 
[BB06] characterization, the corresponding classical Meixner classes are (i) Gaussian (ii) Poisson 
(iii) binomial (iv) gamma (v) hyperbolic secant (vi) negative binomial. The only difference is that 
in the binomial case, — - c e N, but for the free binomial one can take any — 1 < c < 0. See also 
Remark 8. 

Another description of the free Meixner distributions [BB06, Ans07], again with a slightly different 
normalization, is in terms of the i?-transform: for the distribution with parameters b, c, mean m and 
variance t, 

(1) R ^ - m = t + b(R(u) -m) + (c/tMRiu) - mf. 

u 

Note that the second paper quoted above uses the combinatorial /^-transform 1Z(u) = uR(u). 
Remark 5. Our main object it the operator 

LvW = (I ® n)[df\ = [ f{x) ~ fiy) d^y). 

Jr x — y 

We will consider three versions of it. First, L M as an operator on the vector space of polynomials, 

n-l 

Lfj, : R[x] -> R[x], L^] = J2 m n-k-ix k . 

k=0 

Like other deformations of the derivative operator, it lowers the degree by I. See also Lemma 7. 

Second, L M is a linear operator on the vector space of continuously differentiable functions, map- 
ping it to continuous functions, 

: C (1 \R) -> C(R). 

On this space, it has eigenfunctions 



z — x 



z — X 



— = G»(z) — . 
z — X z — X 



For real z, the eigenvalue is irH[fi] (z), where H is the Hilbert transform, whenever this quantity is 
finite. 

The third version is described in the following proposition. See also Lemma 5. 

Proposition 1. Let fi be a compactly supported measure with density w e L 3 (IR, dx). The oper- 
ator is a bounded operator from L 2 (fi) to and from L°°{p) to L 2 {ji). It is a (possibly) 
unbounded, densely defined, closed linear operator on L 2 (/i). 
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Proof. Denoting \\H\\ LP ^ dx) ^ LP{ ^ dx) = C p , 

p[H[p] 2 ] = \\H\tfw\\ x < \\H\ii]\\l \\w\\ 3 < C\ IMI3 • 

So for / G L 2 (p), 

»[\fm\] < vMwm 2 ]<c 3 \\w\\i /2 /i[/ 2 ] i/2 - 

Also, 

p[\H[fp]\] = IW^MIi < \\H\flAh/2 IMIs < C3/2 U/HI3/2 \\w\\ 3 . 



But 



ii/-ii 3/2 = m 3/2 M 3/2 ; < m 3/2 - 3/4 



2/3 



2/3 



'W 



3/4II 2 / 3 



4/3" 

II r-2 M 1 / 2 || 3M 1 / 6 r f 21 1/2 II 1 1 1/2 

= ||/ W\\ \\W 6 \\ = fi[f ] 1 \\w\\ 3 ' . 

Thus 

m[fti\\<cz,2\\w\\T m/ 2 ] i/2 . 

It follows that all the quantities in the following equation are well defined, 

LM\=*f(x)H\pi\(x)-*H\Jii\(x\ 
and ||^|| L 2 (m) ^ l i (m) < 7r(C 3 + C 3/2 ) |MI 3 /2 - Similarly, for / e 

Mtf[//i] 2 ] 1/2 = H^HII' 2 < ||^[/^] 2 || 3 /2 M\T = IWdlla IMlf < c 3 ll/IL IMlf 

and HMl-^hl^) < 27rC 3 ll^lll /2 - 

Finally, let f n -> / in L 2 (^) and # n = £,„[/„] -> # in L 2 (/i). Then -> L M [/] in L l (p), so 
g = L M [/], and L M : L 2 (yu) — > L 2 (p) is closed. □ 

3. The free Meixner characterization 
Theorem 2. For p with all moments finite, consider a "Sturm-Liouville" -type operator of the form 

Q» = p( x ) l I + q( x ) L », 

where p(x), q(x) are polynomials. It has a polynomial system of eigenfunctions only if p is a free 
Meixner distribution. In that case, the conjugate variable for p is ~^y- 

Proof. Clearly 

Q[x] = q(x). 

Thus to have a polynomial eigenfunction of degree 1, we need q to have degree at most 1. This 
implies that p has degree at most 2. Thus denote 

p(x) = a + bx + ex 2 

and 

q(x) = d + ex. 

Next, we compute 

Q[l] = 0, 
Q[x] — (d + ex) 
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and 

Q[x n ] = ( a + bx + cx 2 )(x n ™ 2 + ...) + (d + ex)(x n - 1 + ...) = (c + e)x n + .... 

Thus the eigenvalues (if any) are for n = 0, e for n = 1, and c + e for n > 2. Therefore the 
polynomial system of eigenfunctions can be taken to be 

{1, x + a, x n + f3 n x + 7„, n > 2} 

for some a, f3 n , 7„. We compute 

<5[x + a] = d + ex = e(x + a), 

thus 

(2) ea = d. 
Similarly, for n > 2 

n— 2 n—k—2 / n—l 

Q[x n + (3 n x + 7„] = (a + 6x + ex 2 ) ^ x fe ^ m i m n _ fc _ i _ 2 + (rf + ex) I ^ x k m n _ k _ 1 + (3 n 

k=0 i=0 \fc=0 

= (c + e)(x n + /3 n x + 7„). 
Comparing coefficients, we get for k = 

n-2 

(3) a ^ mim„_i_ 2 + dm n _i + d/3 n = (c + e)7„, 
for fc = 1 

n— 3 n— 2 

(4) a ]T mim n _i_ 3 + 6 ^ mim n _i_ 2 + dm n _ 2 + em n _i + e/5„ = (c + e)(3 n , 

i=0 i=0 

while for 2 < k < n — 1 

n—k—2 n—k—l n—k 

a ^2 m i m n-k-i-2 + b ^ m i m n-k-i-i + c^2,mim n - k -i + c?m ra _ fc _i + em»_ fc = 0. 

i=0 i=0 i=0 

Here the empty sums are understood to be zero. Changing variables from n — k to n, we see that 

for n > 1 

n— 2 n—l ra 

(5) a ^ m i m n-i~2 + m i m n-i-i + c ^ m,im n ^i + dm„_i + em„ = 0. 

i=0 i=0 i=0 

Note that for n = 0, the corresponding term is c + e. In terms of the moment generating function 
M(z) of fi, equation (5) gives 

az 2 M(z) 2 + bzM(z) 2 + cM(z) 2 + dzM(z) + eM(z) =c + e. 

In terms of the Cauchy transform 



G,(z)= [ —dfi(x) = -M(l/z), 
Jm. z ~ x z 



(treated as a formal power series if necessary) the relation is 

(6) (a + bz + cz 2 )G(z) 2 + (d + ez)G(z) - (c + e) = 0. 



BOCHNER-PEARSON-TYPE CHARACTERIZATION OF THE FREE MEIXNER CLASS 

So 

-{d + ez) - y/(d + ez) 2 + 4(c + e)(a + bz + cz 2 ) 
^ ~ 2(a + 6i + cz 2 ) ' 

This means that on the support of /i 

: (d + ex) 2 + 4(c + e)(a + bx + cx 2 ) < 0} U {x : a + bx + cx 2 = 0} , 

the conjugate variable of \x is 

d + ex 



a + bx + cx 2 

On the other hand, in terms of the i?-transform determined by 

G (r(z) + ^)=z, 

equation (6) states that 

cuR(uf + {bu + (2c + e))R(u) + (au + b + d) = 0, 

or 

(7) - (2c + e)^ - — = a + bR(u) + cR(u) 2 , 

u u 

which is equation (1). 

For future reference, we also compute the eigenf unctions. Denoting 

oo 
n=2 

we get from equation (4) that 

cB(z) = az 3 M(z) 2 + bz 2 M(z) 2 + dz 2 M(z) + ez{M(z) - 1) = (c + e)z - czM(z) 2 - ez, 
so either c = or 

B{z) = z{\ - M(z) 2 ). 

Similarly, for 

oo 
n=2 

we get 

(c + e)C(z) = az 2 M(z) 2 + dz(M(z) - 1) + dB(z) = zM(z)(azM(z) + d(l - M{z))) 

ifc^O. □ 

Proposition 3. If fi is a normalized free Meixner distribution with parameters b and c > — 1, the 
corresponding operator Q from Theorem 2 has polynomial eigenfunctions for the unique (up to a 
factor) choice of 

p(x) = 1 + bx + cx 2 , q(x) = —(b + (1 + 2c) x). 
If c — — \, we additionally require that 6 = 0. 
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Proof. By shifting and re-scaling the measure, we may assume that mi = 0, m 2 = 1. From 
equation (5), 

b + 1m\c + d + m\e = 

and 

a + 2mib + (2m 2 + m\)c + mid + m 2 e = 0. 
So under this normalization, 

b + d = 0, a + 2c + e = 0. 
In equation (7), this corresponds to 

a R(u)_ = a + bR ,^ + cR , u \2 

u 

If a = 0, then R(u)(b + cR(u)) = 0, which corresponds to constant R(u) and \x being a delta 
measure. If a ^ 0, we may assume without loss of generality that a > 0. In this case 

= 1 + (b/a)R(u) + (c/a)i?( M ) 2 , 

u 

which is equation (1) in standard form (for mean zero and variance one). In particular, we know 
that R(u) is an /^-transform of a positive measure for any b and for c/a > —1, in other words for 

a + c > 0. 

Note that this implies the eigenvalue c + e = — (a + c) < 0, while e = — (a + 2c) may be positive, 
negative, or zero. 

For e = — (a + 2c) = 0, equation (2) additionally implies that b = —d = as well. □ 
Example 1. We thus get a correspondence between operators 

(py'Y + qy' 

and 

pLl + qL^ 

having polynomial eigenvalues, given by 

q W 

- = — = 

p w 

Note that the operators have slightly different forms, but this is the appropriate correspondence. 
Note also that the parameter ranges in parameterizations in Remarks 1 and 4 coincide in some but 
not all cases, and that classically, there is no analog of the restriction at the very end of the proof 
above. 

The special classes of measures are: semicircular for 

Q = L 2 ^ — xL^ 
with eigenvalue —1, Marchenko-Pastur for 

Q — (a + bx)^ — (b + ax)L^ 

with eigenvalue —a, Bernoulli for 

Q — (a + bx — ax 2 )^ — {b — ax)L /Ji 
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with eigenvalues a, and for general a + c > 0, 

Q = ( a + bx + cx 2 )Ll - (b + (a + 2c)x)L /M . 

If a + 2c = 0, or in other words (c/a) = —\, then also 6 = from the condition ea = d. The 
corresponding distribution is the arcsine law, with 

Q = (2 - x 2 )Lj 

and eigenvalues 0, — \. More generally, in category (iii) or Remark 4, the only values of a, (3 with 
a + (3 = 2 producing polynomial eigenfunctions are a — (3 — 1. Note that classically, these 
parameter values correspond to the uniform distribution, and Legendre polynomials. 

Remark 6. A very different characterization of the free Meixner class in terms of certain operators 
mapping polynomials to themselves appears in [BL09]. It would be interesting to see if there is a 
relation to our results. Note however that the classical version of the results of that paper involves 
the Meixner and not the Bochner class. 

Lemma 4. The only operators of the type in Theorem 2 with polynomial eigenfunctions orthogonal 
with respect to p correspond to the semicircular distributions. 

Proof. By re-scaling, we may assume that p has mean zero and variance one. If the eigenfunction 
polynomials are orthogonal with respect to p, in particular they are centered with respect to it. So 

p[x + a] = m 1 + a = m\ + (d/e) = 

implies a = 0, and from (2), = d = —b. Also 

fi[x 2 + /3 2 x + 72] = m 2 + mif3 2 + 12 = m 2 - 2(c/ c)m\ + - — ^— ^ = 0, 

c + e 

so a = — (c + e) and c = 0. Thus p is the semicircular distribution. □ 

Remark 7 (Free probability). Out of the five Pearson classes, only the first three correspond to 
orthogonal polynomials and the Bochner scheme. Interestingly, exactly the three classes corre- 
sponding to these under Theorem 2, and not the rest of the free Meixner distributions, have an 
interpretation in free probability theory. The semicircle law is the "free Gaussian", since it ap- 
pears as the limit in the free central limit theorem. Similarly, the Marchenko-Pastur law appears 
as the limit distribution in the free version of the Poisson limit theorem. Finally, the free binomial 
distribution is the sum of freely independent Bernoulli variables. See [VDN92] for more details. 

Remark 8 (Random matrices). Let V be a potential such that 

(8) VeC {2) (R), V is convex, / e v ^dx < 00. 

Jr 

Let X be an Hermitian n x n random matrix [Dei99] distributed according to 

(9) ^e' trV ^dX. 

Then as n — > 00, the spectral distribution of X converges to a compactly supported measure p such 
that 

(10) V = 2ttH[p\. 

So if w(x) dx = e~ v ^ dx, then — — = H, the conjugate variable for p. Of course, this was one of 
the original motivations for Voiculescu's definition of the conjugate variable. 
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As shown above, this is also exactly the correspondence between the Pearson and the free Meixner 
classes, although the relevant ranges of parameters do not coincide in all cases. V with V'(x) = 
a+^cxi satisfies condition (8) only if it is an affine transformation of V = x 2 /2, and if X is in the 
Gaussian unitary ensemble 

Z 

then its spectral distribution converges to the semicircular distribution. But also, for a > 1, — (a — 
1) log a; + x satisfies (8) on (0, oo), and if X is in the Wishart ensemble [Mui82] 

^detX^V^W dX 

Z 

(note that in the n = 1 case, this is the gamma distribution), then its spectral distribution converges 
to the Marchenko-Pastur law with parameter a. Similarly, for a, (3 > 1, (a — 1) logx + ((3 — 
1) \og(x — 1) satisfies (8) on (0, 1), and if X is in the Jacobi ensemble [CC04, Col05] 

\ det(J - X)"' 1 det X^I^xki dX 
Z 

(which is a matrix version of the beta distribution), the spectral distribution converges to a free 
binomial law with parameters a, (3. Note that these parameter values correspond precisely to the 
distributions being absolutely continuous. Finally, the remaining two distributions in the Pearson 
class do not correspond to convex potentials, but still satisfy the hypothesis of Theorem 1.3 of 
[ST97, Chapter I], and so have a unique equilibrium measure, which however is not necessarily 
determined by equation (10). 

For a,/3 < 1, one also has matrix models, but they are less canonical. Namely, let P be an n x n 
projection matrix of rank k < n, and X be a GUE matrix as above. Then the matrix 

W = XPX 

has the Wishart distribution with parameters n, k, and as n — > oo while | — > a, the spectral 
distribution of W converges to the Marchenko-Pastur distribution with parameter a, < a < 1 
[Mui82]. Moreover, if W 1: W 2 are Wishart matrices with parameters n, k x and n, k 2 , respectively, 
with ki,k 2 < n, k\ + k 2 > n, then 

(Wi + w 2 )~ 1/2 w 1 (w l + w 2 y 1/2 

is well defined, and as n — > oo, ^ — > a, — > (3, its spectral distribution converges to the free 
binomial distribution with parameters < a, f3 < 1, a + (3 > 1 [CC04]. There is also another 
model for this distribution involving a product of random projections [Col05]. 

While there exist abstract characterizations of, for example, the Wishart distribution [BW02], we 
are not aware of a Bochner type characterization for these matrix distributions for 2 < n < oo. 

We close the section with the following analog of the operator i(y' + being self adjoint with 
respect to w(x) dx; note the absence of boundary conditions. 

Lemma 5. For fi with density in L 3 (R, dx), the operator A = i(L^ — ^H\jj]) is a self-adjoint 
operator on L 2 {ji) with dense domain V — {/ e L 2 {ji) : A[f] G L 2 (/i)}. 
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Proof. By Proposition 1, the operator A maps L°°(fx) to L 2 (fx), so the domain V is dense. Suppose 
that there are g,h £ L 2 (fx) such that for all / G V, 

fx[A[f]g] = fi[fh]. 

For all / G L°°(fx), 

fx[L[f}g] = (fx ® fx)[(g ® I)df] = (fx ® //)[%/) - (/ ® /)£><?] 
= /U [^/]-(/i®/i)[(J® /)£><?] 

= - (// ® //)[(/ ® /)<%] = fx[Hgf] - fx[fL[g}}. 

This implies that the operator A is symmetric and, for such /, 

fx[fh]=fx[A[f]g]=fx[fA[g]\. 

Here we have also used the fact that A maps L 2 (fx) to L}(fx). It follows that h = A[g] in so 
G £ 2 (/x) and g EV. Therefore A is self-adjoint. □ 

4. Extensions 

4.1. Higher-order operators. The operator 

fc=0 

with 

Pfc(x) = a k x k + ... 
has polynomial eigenfunctions if and only if 

n n 
k=0 k=0 

However, unlike in the second order case, where the complete description of free Meixner distribu- 
tions is available, it is not clear which solutions of this equation are Cauchy transforms of positive 
measures. Such analysis is related to the study of free exponential families, see [Bry09]. 

4.2. The two measure case. 

Remark 9. In [Ans09a], with motivation from the two- state free probability theory (also called the 
c-free theory), we introduced the c-free conjugate variable H^y. for a pair of measures fx, v, 
(which may not exist) is determined by 

(fx®v)[df]=fx[H^f\. 

Explicitly, under appropriate conditions on fx, v, 

note however that H^ u may be defined even if H^, H u are not. 
Lemma 6. For any fx, u for which H^ u is well defined, the operator 

is symmetric with respect to the inner product induced by fx. 
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Proof. 



^[LM-H^L^g] = (p®v)[(dL,[f])(g® I)} - ^[H^L^g] 

= -(p ® u)[(I ® L^[f])dg] = -u\L^f]LM- 



□ 



Remark 10. The map p = $p^[v] is the left shift of the Jacobi parameter sequences of the mea- 
sure: it takes a measure with Jacobi parameters {(/3 , Pi, P21 ■ ■ ■ , ), (71, 72, 73, • • •)} to tne measure 
with Jacobi parameters {(/?, (3 , /5i , /? 2 ,...,), (7, 71, 72, 73, •••)}> and so is an inverse of "coefficient 
stripping" [DKS09]. Equivalently, 

M~\z) = l-(3z- 1 z 2 M l/ (z) 

or 



z-P- 1 G v {z)' 

The map $ = $ 01 was defined in [BN08]. In particular, the free Meixner distributions are exactly 
the image under <3> of general semicircular distributions. One can express $/3, 7 [H through <3>[z/] 
using the operation of boolean convolution, see [Ans09b]. 

Lemma 7. A monic polynomial family {A n } such that 

Li/[A n \ = A n _i 

is orthogonal with respect to some p if and only if v is a semicircular distribution and p = ^^[v] 
is a free Meixner distribution. 

Proof. For any p, v, there is a unique monic polynomial family such that L u [A n ] = A n _i and 
p[A n ] = for n > 1. These are the c-free Appell polynomials A%" investigated in [Ans09a], 
where the statement about their orthogonality was also proved. □ 

Lemma 8. p — $ [u] if and only ifH^ u = x. More precisely, for compactly supported p, v, for any 
continuously differentiate f, 

(is®p)[df} = [xf\ 

if and only if p = More generally, H^ v = ^(x — (5) if and only if p = <J>g j7 [i/]. 

Proof. This result was proved in [Ans09a] for polynomial /. Here is a more direct analytic proof. 
If p — <&[i/], then 

G„(z) = z 



Therefore 



(ID 



d- 



z — X 



V 



z — X 



p 



z — X 



z — 



= zG^z) 



p 



X 



G,{z) 
- p[l] = p 



G,{z) 
1 



x- 



X 



Since p, v are compactly supported, this implies that for any n, {y <8> p) [dx n ] = p [x n+1 ]. By 
the Weierstrass theorem, M[x] is dense in C^supp (//)] with the uniform norm on /'. Note that if 
P' n f uniformly on a compact set, by adjusting the constant term we can arrange to have P n — > / 
uniformly. Now take / continuously differentiable. The functional / 1— > (u <g> p) [df] is continuous 
with respect to the uniform norm on /', and the functional /•—>// [xf] is continuous with respect to 
the uniform norm on /. Therefore they coincide on any continuously differentiable function. The 
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converse follows by running equation (11) backwards, and the proof of the more general statement 
is similar. □ 

Proposition 9. For polynomial p, q, the operator 

Q = p(x)L u L lJL + q(x)L fM , 
has a polynomial system of eigenfunctions if and only if 

(az 2 + bz + c)M^(z)M u (z) + (dz + e)M^(z) = c + e, 

where p{x) = a + bx + cx 2 , q(x) = d + ex, and M^(z), M v [z) are the moment generating functions 
of [i, v, respectively. 

The proof is similar to that of Theorem 2. 

Remark 11. The conclusion of Proposition 9 is equivalent to 

(c + e)W i '"(u) = cK v (u) -du- au 2 - buR v [u) - c1Z v {u) 2 , 

where W is the combinatorial i?-transform of v, and W 1 ^ is the i?-transform in the two-state free 
probability theory [Ans09a]. 

Example 2. If b = c = 0, and taking without loss of generality e = —1, we get 

M^z)' 1 = l- dz- az 2 M u (z)i 

in other words \x = $<2,aM- Recall that by Lemma 8 these are exactly the cases when H^ u is a 
polynomial of degree one. Thus for arbitrary v, one can choose a positive ji so that Q has polyno- 
mial eigenfunctions. Moreover, in this case the eigenfunctions are orthogonal. In particular, this is 
the case for v a shifted semicircular distribution and /i the corresponding free Meixner distribution. 

Proposition 10. For p, = Qp^v], the monic orthogonal polynomials P% with respect to \x are 
eigenfunctions for the operator 

Q = L l/ L II - 1 ~ 1 (x-{3)L fl 
with eigenvalue — r y~ 1 forn > 0. A posteriori, this operator is simply 

Q[/] = 7- 1 (-/ + M/l). 

Proof. It is a classical result about monic orthogonal polynomials of the first and second kind that 

^ra = (/®^)[9p^]=^ 1 . 

So using Lemma 8, 

— _A II P" II 2 

— <Jnk H^n-l Hi, • 

But using notation from Remark 10, 

rai = ryi---7»-i = 7||^i|i;, 

so we conclude that 

MOMW]] - h^l,[pz\)p£\ = -r l K k \\pX 

and 

forn>l. □ 



16 M. ANSHELEVICH 

The following result should be well-known, but we have not found it in the literature. The particular 
cases involving Chebyshev polynomials of the first kind, and the free Meixner polynomials [Ans03, 
Theorem 4], are indeed well known. It follows from the preceding proposition, but a direct proof is 
also straightforward. 

Corollary 11. Let fibe a measure with Jacobi parameters {(ft, f3 , fli, fi 2 , • • • , ), (7, 7i, 72, 73, •••)}> 
v the "once-stripped" measure with Jacobi parameters {(/?o, Pi, #2, ■ ■ ■ , ), (71, 72, 73, • • •)}> an d r 
the twice-stripped measure with Jacobi parameters f3 2 , ...,), (72, 73, • • •)}• Then the corre- 
sponding monic orthogonal polynomials satisfy 

xP^x) = PX +1 (x) + (3P^(x) + iP^x). 
Example 3. If v is a Bernoulli distribution, 

v = aiSfr + a 2 5p 2 

with ct\ + a 2 = 1, then 

M (z) - ai + a2 - 1 _ ( ai P 2 + a2 ^ l ) z 



l-(3 lZ 1-foz 1 - (A + fo)z + fcfoz 2 ' 
Now take e = —1, d = ai(3 2 + a 2 (3i, c = 1, b = —(Pi + f3 2 ), a = (3i(3 2 , so that 

(az 2 + bz + c)M v (z) + (dz + e) = = c + e. 
It follows from Proposition 9 that for Bernoulli v, for any \x there is an operator of the form 

Q = (a + bx + x 2 )L u L ll — (x — d)L lx 
with polynomial eigenfunctions. 

Proposition 12. Among all operators pL u L^ + qL^ with polynomial eigenfunctions, only the ones 
in Proposition 10 have orthogonal eigenfunctions. 

Proof By re-scaling, we may assume that p has mean zero and variance one. Denoting by {m v k } 
the moments of v, it follows from Proposition 9 that 

(12) b + cm u 1 +d = 
and 

(13) a + bm\ + c + cm^ + e = 0. 
Also, using the notation from the proof of Theorem 2, 

cB(z) = cz(l - Mp{z)M v (z)) 

and 

(c + e)C(z) = zM^z)(azM u (z) + d(l - M v (z))). 

Thus 

p[x + a] = + a = + (d/e) = 
implies a = 0, so d = and b = —cm". Also 

fi[x 2 + (3 2 x + 72] =«i2+ m iP2 + 72 = m 2 - (c/c)m^(m^ + mj) + - — ^— ^ = 0, 

c + e 

so a = dm^ — (c + e) = — (c + e). It follows from equation (13) that bm\ + cm v 2 = = 
c (mj — {yn\) 2 ), so unless z/ is a delta measure, c = and 6 = 0. 
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If v — 5$ with M v (z) = Tzji> then using Proposition 9 and the calculations above, 

'— (c + e)z 2 — c£z + c 



and 



(c + e)(l-£s) 



f e ) Mn(z) = c + e 



— (c + e)z 2 — c£z + c + e — e£z 1 — £z — z 2 
It follows that pi is a (Bernoulli) free Meixner distribution, = pL^-i = $ [^] = 



□ 
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